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Particle acceleration in collisionless shocks is believed to be responsible for the production of 
cosmic-rays over a wide range of energies, from few GeV to > 10 20 eV, as well as for the non- 
thermal emission of radiation from a wide variety of high energy astrophysical sources. A theory 
of collisionless shocks based on first principles does not, however, exist. Observations of 7-ray 
burst (GRB) "afterglows" provide a unique opportunity for diagnosing the physics of relativistic 
collisionless shocks. Most GRBs are believed to be associated with explosions of massive stars. Their 
"afterglows," delayed low energy emission following the prompt burst of 7-rays, are well accounted 
for by a model in which afterglow radiation is due to synchrotron emission of electrons accelerated in 
relativistic collisionless shock waves driven by the explosion into the surrounding plasma. Within the 
framework of this model, some striking characteristics of collisionless relativistic shocks are implied. 
These include the generation of downstream magnetic fields with energy density exceeding that of the 
upstream field by ~ 8 orders of magnitude, the survival of this strong field at distances ~ 10 10 skin- 
depths downstream of the shock, and the acceleration of particles to a power-law energy spectrum, 
dlogn/dloge ~ —2, possibly extending to 10 20 eV. I review in this talk the phenomenological 
considerations, based on which these characteristics are inferred, and the challenges posed to our 
current models of particle acceleration and magnetic field generation in collisionless shocks. Some 
recent theoretical results derived based on the assumption of a self-similar shock structure are briefly 
discussed. 



I. INTRODUCTION 



Due to the low densities characteristic of a wide range of astrophysical environments, shocks observed in many 
astrophysical systems are collisionless, i.e. mediated by collective plasma instabilities rather than by particle-particle 
collisions. For example, collisionless shocksplayan important role in supernova remnants [e.g.Q, jets of radio galaxies 
[e.g. [3, gamma-ray bursts [GRBs, e.g. 0, Sill, and the formation of the large scale structure of the Universe [e.g. 
01 SHI- Although collisionless shocks have been studied for several decades, theoretically and experimentally, in space 
and in the laboratory, a self-consistent theory of collisionless shocks based on first principles has not yet emerged [see, 
e.g., comments in|lOj. 

Observations of GRB "afterglows," the delayed low energy emission following the prompt 7-ray emission, provide 
a unique probe of the physics of collisionless shocks. Current understanding suggests that the afterglow radiation 
observed is the synchrotron emission of energetic, non-thermal electrons in the downstream of a strong collisionless 
shock driven into the surrounding interstellar medium (ISM) or stellar wind. These collisionless shocks start out highly 
relativistic, with shock Lorentz factor T ~ 100 on time scale of minutes after the GRB, and gradually decelerate to 
r ~ 10 on a day time scale and r ~ 1 on a month time scale. This allows one to probe the physics of the shocks over 
a wide range of Lorentz factors. 

In this talk, I focus on what we have learned about collisionless shocks from the observations of GRB afterglows, 
and on the challenges to our theoretical understanding of the relevant physics. A brief introduction to GRBs and their 
afterglows is given in § [D] The dynamics of shock expansion during the afterglow phase is described in § III Al and 
the model for emission of afterglow radiation is discussed in § III Bl The collisionless nature of the shock is discussed 
in S liTTl 

The theoretical challenges are discussed in § IIIII Our current understanding of, and open questions related to, the 
generation of magnetic fields are discussed in § IIII Al where it is pointed out that the main challenge is the survival 
of the field at distances much larger than the skin depth downstream of the shock. Our current understanding of, 
and open questions related to, particle acceleration are discussed in § IIII Bl One of the open questions related to 
particle acceleration is whether GRBs may produce the observed ultra-high energy (UHE), > 10 19 eV, cosmic-rays. 
The challenges posed to all models of particle acceleration by the observations of UHE cosmic-rays, and the arguments 
suggesting GRBs may be the sources of these particles, are discussed in § IIII CI 

In the past few years, a significant effort was invested in 3D numerical studies of collisionless shocks. These studies 
are briefly discussed in § I1VI with an emphasis on what we have learned from the simulations, and on whether 
such simulation are likely to provide answers to the main open questions (for a more detailed summary of numerical 
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simulation results see HH). Recent theoretical results derived based on the assumption of a self-similar shock structure 
are briefly discussed in § [V] A brief summary is given in § I VII 

It should be pointed out that although the discussion in this lecture is motivated by GRB afterglow observations, 
it may be relevant also for non relativistic collisionless shocks, such as shocks in young supernova remnants (SNRs) 
and in the intergalactic medium. In the past few years, high resolution X-ray observations have provided indirect 
evidence for the presence of strong magnetic fields, > 100/iG, in the non-relativistic (v/c ~ 10~ 2 ) shocks of young 
SNRs (e.g. ^3 an d references therein). These fields extend to distances D > 10 17 cm ~ 10 10 skin depths downstream, 
and possibly even > 10 16 cm upstream, of the shock. In resemblance to GRBs, such strong magnetic fields cannot 
result from the shock compression of a typical interstellar medium (ISM) magnetic field, £?ism ~ few fj,G. In SNRs the 
discrepancy is somewhat less severe, and the possibility that these magnetic fields are related to the large scale ISM 
fields cannot be ruled out. If the ISM magnetic fields can be neglected for SNR shocks, then much of the discussion 
presented below applies to such shocks as well. 



II. GAMMA-RAY BURSTS AND THEIR AFTERGLOW 
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Gamma-ray bursts (GRBs) are short, typically tens of seconds long, flashes of gamma-rays, carrying most 

of their energy in > 1 MeV photons. The detection in the 
past few years of "afterglows," delayed X-ray, optical and ra- 
dio emission from GRB sources, proved that the sources lie at 
cosmological distances, and provided strong support for the 
scenario of GRB production described in fig. ^ The energy 
source is believed to be rapid mass accretion onto a newly 
formed solar-mass black hole (or, possibly, neutron star). Re- 
cent observations suggest that the formation of the central 
compact object is associated with type Ib/c supernovae. 

The energy release drives an ultra-relativistic plasma out- 
flow, with Lorentz factor T ~ 10 2 5 . The emission of 7-rays 
is assumed to be due to internal collisionless shocks within 
the relativistic wind, the "fireball," which occur at a large 
distance from the central black-hole due to variability in the 
wind emitted from the central "engine." It is commonly as- 
sumed that electrons are accelerated to high energy within 
the collisionless shocks, and that synchrotron emission from 
these shock accelerated electrons produces the observed 7- 
rays. At still larger distances, the wind impacts on the sur- 
rounding medium. Here too, a collisionless shock driven into 
the ambient gas is believed to accelerate electrons to high energy, leading to synchrotron emission which accounts for 
the "afterglow." Afterglows are much better understood than the prompt 7-ray emission. 
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FIG. 1: The fireball scenario of GRB production 



A. Afterglow I: Dynamics 
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As the fireball expands, it drives a relativistic shock (blast wave) into the surrounding gas, e.g. into the 

interstellar medium (ISM) gas of the galaxy within which the ex- 
plosion occurs. At early times, the fireball is little affected by the 
interaction with the ISM. When the mass of the shocked ISM plasma 
becomes significant, typically exceeding E/Tc 2 where E is the ex- 
plosion energy, the fireball begins to decelerate. At sufficiently late 
time, when the shock radius becomes much larger than the radius of 
the onset of deceleration, most of the fireball energy is transferred 
to the ISM, and the flow approaches a self-similar behavior [l3j . 

The dynamics of the shock can be easily understood by consider- 
ing the flow in the shock frame, where the shock is stationary and 
the upstream plasma approaches it with high speed, corresponding 
to a Lorentz factor T. At the shock, protons are scattered and their 
momentum distribution is isotropizcd. The bulk velocity of the fluid 
is thus reduced, and a significant fraction of the incoming kinetic 
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FIG. 2: The afterglow shock (shock frame). 
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energy is converted to thermal energy. The characteristic Lorentz 
factor associated with the resulting thermal (i.e. random) motion of 
protons is ~ T. In the observer frame, where the shock moves forward with Lorentz factor T, the energy of each proton 
is ~ T 2 m p c 2 . The evolution of T with shock radius R is determined by energy conservation, E ~ T 2 (R)m p c 2 nR 3 
where n is the ambient medium number density. A more detailed analysis gives |l3j 



1/2 / p ^ 1/2 / R N -3/2 



^ ^ \ lQiTnm p c 2 / (no ) 



10 17 cm 



(1) 



where E = 10 53 -Es3 erg and n = lno cm 3 is the number density. £"53 ~ 1 for typical GRBs, and no ~ 1 for typical 
ISM. 

At the ISM rest frame, the shock reaches a radius R at time R/c. The time measured in the shock frame is shorter 
by a factor T, R/Tc. Since the plasma flows away from the shock (in the downstream region) at velocity ~ c/3 ^ij, the 
shocked plasma is concentrated in a narrow shell of thickness R/T (which becomes R/T 2 in the ISM rest frame due to 
Lorentz contraction), with proper density Tn (corresponding to a density T 2 n in the ISM frame). Finally, the radiation 
emitted by the shocked plasma at radius R is delayed compared to that emitted R — by R/c — R/v = R/2T 2 c, here 
v ~ (1 — l/2r 2 )c is the shock velocity (to leading order in 1/F 2 ). A more detailed calculation gives [w| the following 
relation between the radius and the time measured by a distant observer: t bs. (R) = R/4T(R) 2 c, which implies 



S(W ) = 3.2xlO"(^) ,/1 ( i ^) 1/4 c nl . (2) 

Eqs. J3> and @ provide a complete description of the dynamics of the afterglow shock. The self-similar model 
for shock expansion du ring afterglow emission has been tested by direct measurements of the fireball size, through 
diffractive scintillation 16] and through very large baseline interferometry 17], and by indirect size measurements 
using late time radio observations |1S|. The simple self-similar model predictions are in excellent agreement with 
observations. 



B. Afterglow II: Radiation 



The agreement between the dynamics predicted by the blast wave model and the direct measurements of the fireball 
size strongly argue for the validity of this model's dynamics. As explained in § III CI the shock wave is most likely 
collisionless, i.e. mediated by plasma instabilities. The electromagnetic instabilities mediating the afterglow shock are 
expected to generate magnetic fields, and the rarity of binary particle collisions is expected to lead to the generation 
of a non-thermal particle distribution downstream [see, e.g. IjJ, for a discussion of particle acceleration]. Afterglow 
radiation was therefore predicted to result from synchrotron emission of shock accelerated electrons [l9J . The observed 
spectrum of afterglow radiation is indeed remarkably consistent with synchrotron emission of electrons accelerated to 
a power-law distribution (see below) , providing strong support to the validity of the standard afterglow model based 
on synchrotron emission of shock accelerated electrons 

In order to determine the luminosity and spectrum of synchrotron radiation, the strength of the magnetic field and 
the energy distribution of the electrons must be determined. Due to the lack of a first principles theory of collisionless 
shocks, we present in this section a purely phenomenological approach to the model of afterglow radiation emission. 
That is, we do not discuss the processes responsible for particle acceleration and magnetic field generation. Rather, 
we simply assume that a fraction es of the post-shock thermal energy density is carried by the magnetic field, that a 
fraction e e is carried by electrons, and that the energy distribution of the electrons is a power-law, d\ogn e /dloge — p 
(above some minimum energy £0 which is determined by e e and p). e e and p are treated as free parameters, 
to be determined by observations. It is important to clarify here that the constraints implied on these parameters 
by observations are independent of any assumptions regarding the nature of the afterglow shock and the processes 
responsible for particle acceleration or magnetic field generation. Any model proposed for the shock should satisfy 
these observational constraints. 

The parameters e^, e e andp, together with the parameters E and n which determine the shock dynamics, completely 
determine the magnetic field strength and electron distribution (including their temporal and spatial dependence). 
Thus, they completely determine the luminosity and spectrum of synchrotron emission, including its time dependence. 
The power-law dependence of T on time, combined with the power-law dependence of the number of electrons on en- 
ergy, imply that the time dependence of the synchrotron emission (at a given frequency) and its frequency dependence 
(at a given time) are both power-laws. The indices of these power-laws are both determined by p alone. One of the 
most impressive successes of the model described here for afterglow emission is that the observed frequency and time 
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dependence of the afterglow flux both follow power-laws, and that the temporal and frequency power-law indices are 
both consistent with the same value of p. Moreover, the inferred value of p, p w 2, is consistent with the theoretical 
expectations for particle acceleration in collisionless shocks (this is discussed in detail in § IIII B(l . 

In addition to p, the afterglow model is determined by 4 parameters, {E, n, cb, e e }- In principle, all 4 parameters 
may be determined by observations, which provide 4 observables. The synchrotron spectrum can be described as a 
combination of 4 power-laws, with 3 break frequencies: the frequency v a below which the synchrotron self-absorption 
optical depth exceeds unity, the frequency v m where the flux peaks (corresponding to the synchrotron frequency of 
electrons with energy £q), and the frequency v c corresponding to synchrotron frequency of electrons at energy e c for 
which the synchrotron cooling time is comparable to the dynamical time (the time for significant shock expansion). 
Typically, v a < v m < u c , and the specific intensity /„ is a broken power-law, /„ cx v a with a = 2, 1/3, — (p — l)/2, — p/2 
from low to high frequency. Afterglow observations therefore provide 4 observables, v a , v m , v c and the normalization 
of the flux, e.g. j v {y = v m ). This, in turn, allows one to determine the 4 model parameters, {E, n, es, e e }. 

As noted above, the values which are typically obtained for E and n are E ~ 10 53 erg and n ~ 1cm -3 . The inferred 
energy is consistent with that expected based on 7-ray observations, and the inferred density is consistent with that 
typically expected in the ISM. Here we are more interested in the inferred values of es, £e- It is natural to hope that 
the values of eb , e e (and p) are universal since they are determined by the microphysics of the collisionless shock. The 
constancy of p and of e e among different bursts is strongly supported by observations. Universal values of p and e e , 
p w 2 and e e « 0.1, typically inferred from most optical afterglows, are also inferred from the clustering of explosion 
energies [2(j and from X-ray afterglow luminosity [2ll 1221 ] . The value of is less well constrained by observations. 
However, in cases where eb can be reliably constrained by multi waveband spectra, values close to equipartition, 
€b = 0.01 to 0.1, are inferred [e.g.flj|. Observations are consistent with the values of eb, e e and p being independent 
ofT. 

It should be pointed out here that Eichler & Waxman |2^| have shown that observations determine e e and eb (and 
also E and n) only up to a factor /, the fraction of electrons accelerated, where m e /m p < f < 1. However, it is 
expected that / is not very small, / > 1/10 pflj . 

C. Why collisionless shock? 

As explained in the preceding sub-sections, afterglow observations strongly support a model where the afterglow is 
produced by a relativistic shock wave driven by the explosion into the surrounding medium, and where the radiation 
is due to synchrotron emission of high energy electrons, following a power-law energy distribution, radiating in the 
post-shock magnetic field. These conclusions are independent of assumptions related to the nature of the shock and 
to the processes responsible for particle acceleration and magnetic field generation. The nature of the shock is the 
main subject discussed in this sub-section. 

The scattering of particles at the shock, leading to isotropization of momenta in the downstream, can not be 
mediated by binary particle collisions. Let us first consider shock mediation through binary Coulomb collisions. At the 
shock frame, the number density of particles is Tn (both in the up- and downstream), and the cross section for Coulomb 
collisions may be estimated as ttcP with d ~ e 2 /Tm p c 2 . The resulting mean-free-path, A ~ 1/Tnird 2 ~ 10 31 TnQ 1 cm, 
is many orders of magnitude larger than the size of the system, ~ 10 17 cm (in fact, it is larger than the size of the 
observable universe, ~ 10 28 cm). Coulomb collisions can not therefore mediate the shock. The mean free path for 
nuclear collisions, A ~ 1/Tna pp ~ 10 25 (rrio) cm, is also many orders of magnitude larger than the size of the 
system. This implies that nuclear collisions can not mediate the shock either. Note, that we are considering here the 
scattering of protons (rather than electrons), since they carry most of the momentum of the incoming plasma. 

The scattering of particles at the shock is most likely mediated by collective plasma instabilities (involving a 
macroscopic number of particles). As the fast, expanding downstream plasma tries to propagate through the upstream 
plasma, which is at rest, electromagnetic instabilities develop, generating electric and magnetic fields in the plasma, 
which deflect the particles and tend to isotropize their momentum distribution. These instabilities develop on a time 
scale comparable to the inverse of the (relativistic) plasma frequency, t ~ l/u> p . Since the particle distribution is 
expected to be isotropized on a time scale of 1/uj p , the width of the shock (the size of the region over which particles 
are being scattered) is expected to be of the order of several skin depths, few times c/oj p , where 

c/uj p = c(4TrTne 2 /Tm p y 1/2 = c{Ai:ne 2 /m p )~ 1/2 ~ 10 7 n 1/2 cm. (3) 

Note, that in the shock frame both the number density and the "effective" mass of protons are larger by a factor T 
than their values in the ISM, which implies that the plasma frequency is similar to that of the ISM. 

A shock which is mediated by collective plasma instabilities, instead of by binary particle collisions, is termed 
"collisionless." Afterglow shocks, as shocks in many other astrophysical systems, are collisionless due to the low 
densities characteristic of these systems. 
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The following comment is in place at this point. The arguments presented above indicate that a collisionless 
shock would be formed, with a width of the order of several skin depths. This conclusion is supported by numerical 
simulations (see § IIVI for discussion) . It should however be kept in mind that a self-consistent theory of collisionless 
shocks does not yet exist, and that the nature of the shock may thus be different than described above. In this 
context it is important to clarify that the conclusions inferred in § III Al regarding the existence of a blast wave and its 
dynamics, and in S III Bl regarding the strength of the post shock magnetic field and the electron energy distribution, 
are inferred directly from afterglow observations and are independent of any assumptions regarding the nature of the 
shock and the processes responsible for particle acceleration or magnetic field generation. One may not rule out, of 
course, the possibility of explaining afterglow observations with a different model, where a shock wave does not form 
or where radiation is not produced by synchrotron emission. Such an alternative models have not, however, been put 
forward yet. 

III. CHALLENGES 
A. Magnetic field generation 

Afterglow shocks are highly "non-magnetized:" The ratio of magnetic field to kinetic energy flux ahead of the shock 
is very small, t/B iUp /nm p c 2 ~ 10 -10 , where Us, U p is the magnetic energy density in the upstream, and B up ~ 3[i G 
is typical to the ISM. This stro ngly suggests that the shock structure is determined by the upstream density and the 
shock Lorentz factor alone [e.g.|27|. It is reasonable to assumed that the upstream magnetic field does not play a role 
in the determination of the shock structure, since the ratio of the (relativistic) cyclotron frequency of the thermal 
protons, ul ~ eB /Tm p c, to the plasma frequency is (o>l/w p ) 2 ~ UB,up/nm p c 2 ~ 10~ 10 , which implies that thermal 
protons in the shock frame are not affected by the compressed upstream magnetic field on a the dynamical time 

The downstream magnetic field, implied by afterglow observations, is close to cquipartition, i.e. it's energy density 
is higher than that of the upstream field by a factor of ~ 10 8 . Such near equipartition magnetic field may conceivably 
be produced in the collisionless shock driven by the GRB explosion by electromagnetic (e.g. Weibel) instabilities [e.g. 
HI 124. l25l l26| . The main challenge associated with the downstream magnetic field is related to the fact that in order to 
account for the observed radiation as synchrotron emission from accelerated electrons, the field amplitude must remain 
close to equipartition deep into the downstream, over distances ~ 10 10 c/ujp: At t ~ 1 day the magnetic field must 

be strong throughout the (proper) width ~ 2Tct ~ R/T ~ 10 17 cm while c/uj p ~ 10 7 n cm. This is a challenge 
since electromagnetic instabilities are believed to generate (near-equipartition) magnetic fields with coherence length 
L ~ c/Up, and a field varying on such scale is expected to decay within a few skin-depths downstream j23|. This 
suggests that the correlation length of the magnetic field far downstream must be much larger than the skin depth, 
L 3> c/iOp, perhaps even of the order of the distance from the shock 0jI13- 

B. Particle acceleration 

The mechanism which is commonly believed to be responsible for the production of non-thermal distributions of 
high energy particles in many astronomical systems (e.g. planetary bow shocks within the solar wind, supernovae 
remnant shocks, jets of radio galaxies, GRB's and possibly shocks involved in the formation of the large scale structure 
of the universe) is the diffusive (Fermi) acceleration of charged particles in collisionless shocks. In this process, high 
energy particles are scattered back and forth across the shock discontinuity. Since the flow at the shock is converging 
(that is, the fluid on each side of the shock discontinuity "sees" the fluid at the other side as approaching it), each 
time the particle is scattered from one side of the shock to the other it gains energy. Since the scattering is mediated 
by interactions with plasma waves, incorporating a macroscopic number of particles, the scattered particle can reach 
very high energies. One may consider an analogy where a light ping-pong ball is scattered back and forth between 
approaching trains, gradually gaining energy with each collision. The accelerated particle gains energy until it escapes 
the shock at the far downstream. 

Despite decades of research, this mechanism is still not understood from first principles (see, e.g., [28, 29] for discus- 
sions of alternative shock acceleration processes) . Particle scattering in collisionless shocks is due to electro- magnetic 
waves. No present analysis self-consistently calculates the generation of these waves, the scattering of particles and 
their acceleration. Most analyses consider, instead, the evolution of the particle distribution adopting some Ansatz 
for the particle scattering mechanism (e.g. diffusion in pitch angle), and the "test particle" approximation, where 
modifications of shock properties due to the high energy particle distribution are neglected [see, however |3(J| . 
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This phenomenological approach proved successful in accounting for non-thermal particle distributions inferred 
from observations. The theory of diffusive particle acceleration in non-relativistic shocks was first developed in 
HI El El Hi, an d was shown to predict a power-law distribution of particle momenta, dn/de cx e~ p , with 



P 



3/3, 



I3 u - fa 



(4) 




FIG. 3: Accelerated particle spectral index, s — 
p + 2, as a function of Fv/c (where T and v are 
the shock Lorentz factor and speed) for various 
equations of state of the fluid. Solid lines show 
the analytic result, eq. ©, and symbols show the 
results of numerical simulations. 



Here (3 U (Pd) is the upstream (downstream) fluid velocity normalized 
to the speed of light. For strong shocks in an ideal gas of adiabatic 
index 7 = 5/3, this implies p — 2, in agreement with observations 
of non-relativistic shocks. 

Observations of GRB afterglows lead to the conclusion that the 
highly relativistic collisionlcss afterglow shocks produce a power- law 
distribution of high energy particles with p = 2.2 ± 0.2 |2ll l22l l35| . 
This triggered a numerical investigation of particle acceleration in 
such shocks [3(J . The index p was calculated under the " test parti- 
cle" approximation for a wide range of Lorentz factors and various 
equations of state [H El El an d references therein] . In particular, 
p was shown to approach the value 2.2 for large Lorentz factors, in 
accord with GRB observations. These studies have assumed rest 
frame diffusion in pitch angle 36] or in the angle between particle 
velocity and shock normal [371 138| . These two assumptions may not 
be valid, for example if large angle scattering prevails [3^. Both 
yield similar spectra for ultra-relativistic shocks |40| . 

Recently, an analytical study of diffusive particle acceleration in 
relativistic, collisionless shocks yielded a simple relation between the 
spectral index p and the anisotropy of the momentum distribution 
along the shock front [4l(. Based on this relation, a generalization 
of eq. 10} was derived for relativistic shocks, under the assumption 
of isotropic diffusion, 



P Pu~Pd 



(5) 



This result is in agreement with previous numerical determinations of p for all shock speeds and fluid equations of 
state (see fig. |3J). In particular, it yields p = 20/9 = 2.22 in the ultra-relativistic limit. 

The results described above indicate that we have a good understanding of particle acceleration in relativistic 
shocks in the test particle approximation. Several major open questions remain. First, the spectral index p depends 
on the form of the diffusion function (unlike the situation in non relativistic shocks), and a value consistent with 
observations, p ss 2, is obtained for isotropic diffusion only. Second, the fraction of particles which arc accelerated 
can not be determined. Finally, and most importantly, a theory that self-consistently describes the acceleration of 
particles and the generation of electromagnetic waves is still missing. The accelerated particles are estimated to carry 
a considerable part of the energy: electrons alone carry ~ 10% of the energy in GRB afterglow shocks and ~ 5% 
of the energy in SNR shocks |42l and references therein], and at least 10% of the energy in SNR shocks must be 
converted into relativistic protons if these shocks are responsible for Galactic cosmic rays [e.g. El]- This implies that 
the accelerated particles are likely to have an important role in generating and maintaining the inferred magnetic 
fields. This conclusion is supported also by the evidence of strong amplification of the magnetic field in the upstream of 
GRB afterglow shocks |4J] , which is most likely due to the streaming of high energy particles ahead of the shock. Since 
the high energy particles are likely to play an important role in the generation of the fields, a theory of collisionless 
shocks must provide a self-consistent description of particle acceleration, which depends on the scattering of these 
particles by magnetic fields, and field generation, which is likely driven by the accelerated particles. 



C. High energy cosmic rays 

The cosmic-ray spectrum extends to energies ~ 10 20 eV |4^, and is likely dominated beyond ~ 10 19 eV by extra- 
Galactic sources of protons. The origin of the ultra high energy (UHE), > 10 19 eV, cosmic rays is a mystery: The 
stringent constraints, which are imposed on the properties of possible UHE cosmic ray sources by the high energies 
observed, rule out almost all source candidates. The essence of the challenge of accelerating particles to > 10 19 eV can 
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be understood using the following simple arguments, which are independent of the details of the acceleration model 

El- 
Consider an astrophysical source driving a flow of magnetized plasma, with characteristic magnetic field strength 
B and velocity v. Imagine now a conducting wire encircling the source at radius i?, as illustrated in fig. The 

~ • L - J il - 1 v - L '- — J — 5 — ^ive of the magnetic flux $ and is therefore given 

■ated by this potential drop would reach energy 
e p ~ (3eBR. The situation is somewhat more complicate in the 
case of a relativistic outflow, where T = (1 — /3 2 ) -1 / 2 ^> 1. In this 
case, the proton is allowed to be accelerated only over a fraction 
of the radius R, comparable to R/T. To see this, one must realize 
that as the plasma expands, its magnetic field decreases, so the time 
available for acceleration corresponds, say, to the time of expansion 
from R to 2R. In the observer frame this time is R/c, while in 
the plasma rest frame it is R/Tc. Thus, a proton moving with the 
magnetized plasma can be accelerated over a transverse distance 
~ R/T. This sets a lower limit to the product of the magnetic field 
and source size, which is required to allow acceleration to e p , 




1 BR 2 
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FIG. 4: Potential drop generated by a non-steady 
outflow of magnetized plasma. 



BR > Te p /e(3. 



(6) 



Eq. [S] also sets a lower limit to the rate L at which energy should be generated by the source. The magnetic field 
carries with it an energy density B 2 /8ir, and the flow therefore carries with it an energy flux > vB 2 /8tt (some energy 
is carried also as plasma kinetic energy), which implies L > vR 2 B 2 . Using eq. |S|we find 
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erg/s. 
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Only two types of sources are known to satisfy this requirement. The brightest steady sources are active galactic 
nuclei (AGN). For them T is typically between 3 and 10, implying L > 10 47 erg/s, which may be satisfied by the 
brightest AGN. The brightest transient sources are GRBs. For these sources T ~ 10 25 implying L > 10 50 5 erg/s, 
which is generally satisfied since the typical observed MeV-photon luminosity of these sources is L 1 ~ 10 52 ergj/s. For 
a more detailed discussion of the arguments suggesting an association of GRBs and UHE cosmic-rays see |42| and 
references therein. 



IV. NUMERICAL SIMULATIONS 



The search for a self-consistent theory of collisionless shocks has led to extensive numerical studies. Particle in cell 
(PIC) simulations were performed in one dimension [e.g. \A8\. in two dimensions [e.er. I49L |H(il[ n in tw o sp atial and three 
velocity dimensions [2D3V; e.g. l27ll5lll52| and in the last few years also in three dimensions |53U54ll55| . In addition to 
shock simulations, homogeneous 3-dimensional simulations of inter-penetrating plasmas have been performed [5ril57l 
IHsj . It is important to note here that the applicability of 1-dimensional and 2-dimensional simulations to the problem 
of collisionless shocks is questionable, since the 1- or 2-dimensional symmetry imposed prevents the evolution of 3- 
dimensional modes, which are likely important for the generation of the shock structure. Similarly, the applicability of 
homogeneous (inter-penetrating) plasma simulations is questionable, since the process of particle acceleration, which 
is likely related to scattering of particles across the shock and which may play an important role in generating large 
scale magnetic fields, may be completely absent in such simulation. 

The 3-dimensional simulations have provided compelling evidence that transverse electromagnetic (Weibel- like) 
instabilities can mediate relativistic collisionless shocks, and that the upstream magnetic field is unimportant if the 
upstream is not highly magnetized. As expected, the width of the shock transition region is found to be a few 
xlO c/ojp [see, e.g.. l55j. These simulations also indicate that near-equipartition magnetic fields are generated in 
the downstream of relativistic shocks propagating in pair (e + e _ ) plasmas, and £b > fn e /m p magnetic fields are 
generated in the downstream of shocks in electron-proton plasmas. Simulations of electron-proton plasmas are forced 
to employ an effective small proton to electron mass ratio, rh p /m e < 20 with present computational resources, and 
the preliminary results thus obtained are not easily extrapolated to more realistic mass ratios. 

3D simulations are limited to very small simulation boxes, and can reliably probe small length scales no larger 
than ~ 100 electron skin-depths, and short time scales no longer than ~ 100 electron plasma times. Obviously, the 
question of field survival and correlation length evolution on length scales 3> c/oj p , which is the main challenge to our 
theoretical understanding, are not yet answered. Similarly, highly energetic particles cannot be contained in the small 
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simulation boxes used, so Fermi-like acceleration processes are suppressed. It appears, therefore, that the numerical 
resources required to answer some of the main open questions by direct 3-dimensional numerical simulations will not 
be available in the near future. 

It is also important to note here, that some published results are based on PIC simulations in stages where the 
boundary conditions strongly modify the plasma evolution. For example, claims that the magnetic fields decay slowly 
or saturate at some finite level remain questionable, until verified by simulations with sufficiently large simulation 
boxes. 



V. SELF-SIMILARITY 



As discussed in § 1111 Al afterglow observations suggest that the characteristic length scale L for variations in the 
magnetic field becomes much larger than the skin depth, L 3> c/uj pi at distances D 3> c/oj p downstream of the shock. 
It is thus reasonable to assume that, at D 3> c/u) p , L is the only relevant length scale, which implies self-similarity 
|ll| . There is no proof that self-similarity will be present whenever the characteristic length scale diverges (or becomes 
infinitesimal). However, the self-similarity assumption is known to be valid for many physical systems in which such 
divergence occurs (see, e.g., [59] for self-similarity in hydrodynamics, and 60] for self-similarity in critical phenomena). 

Self-similarity implies that the plasma configurations at different distances D downstream of the shock, correspond- 
ing to different values of L, are similar to each other. Consider for example the average of the particle distribution 
function over planes perpendicular to the shock normal, < f a > (_D,p), or the magnetic field correlation function, 
Bij (Ax, At, D) =< Bi(x,t)Bj(x + Ax, t + At) > (Here, p is the particle momentum, <> denotes an average over 
planes perpendicular to the shock normal, and we have assumed that such averages depend only on D). Self-similarity 
implies that such averages scale with L, i.e. that 



< f a > (p, L) 



L 



< fa > 



(L/L ) 



-,L>o 



(8) 



B%j (Ax, At, L) = [ A 



2*i 



Aj 



At 



L/Lo' (X/L )* 



(9) 



Here, we have replaced the dependence on D with a dependence on L, as the self-similarity assumption implies 

L oc D. 



Using the Maxwell- Vlasov equations, the similarity indices may be derived [TT| : 

-1 < s B < 0, s t = 1, s p = s B + l, Sf = —&-2s B . 



(10) 



(11) 



These relations imply that the characteristic Larmor radius of energetic particles scales as L, and that the energy 
density of energetic particles in any momentum interval, with the interval scaling as L Sp , scales as the magnetic field 
energy density oc L 2sb . Under the assumption that accelerated particles reach the shock front (and/or are advected 
to the downstream) , the spectrum of accelerated particles must follow 



dn/de(xe- 2/{sB+1) . 



(12) 



The self-similarity assumption suggests that the plasma may be approximately described as a combination of 
two self-similar components: a kinetic component of energetic particles, and an MHD-like component representing 
"thermal" particles. It is likely that the thermal component may be considered as an infinitely conducting fluid, 
in which case sb = and the scalings are completely determined, e.g. dn/de oc e -2 and B oc D°, with possible 
logarithmic corrections. 

The self-similarity assumptions may be tested through their predictions for the evolution of homogenous (time- 
dependent) plasmas, which may be accessible to direct numerical simulations. An inclusion (at the initial conditions) 
of a non isotropic, power-law spectrum of high energy particles, dn/de oc e 2+lp , in homogenous simulations may lead 
to a self-similar evolution in time, with magnetic field evolution following B oc £-( i p+ 4 )/ 2 ('p+ 3 ). 



VI. SUMMARY 



1. GRB afterglows provide a unique probe of relativistic collisionless shocks, with Lorentz factors in the range of 
- 100 to - 1. 
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2. Afterglow shock ares "non-magnetized," in the sense that the ratio of magnetic field to kinetic energy flux ahead 
of the shock is very small, UB, U p/nm p c 2 ~ 10~ 10 , where Us.up is the magnetic energy density in the upstream. 
This suggests that the upstream field plays no role in determining the shock structure. 

3. The downstream magnetic field implied by afterglow observations is close to equipartition, i.e. it's energy 
density is higher than that of the upstream field by a factor of ~ 10 8 . Such near equipartition magnetic field 
may conceivably be produced by electromagnetic (Weibel- like) instabilities. 

4. The main challenge associated with the downstream magnetic field is related to the fact that in order to 
account for the observed radiation as synchrotron emission from accelerated electrons, the field amplitude must 
remain close to equipartition deep into the downstream, over distances ~ 10 10 c/o> p . This is a challenge since 
electromagnetic instabilities are believed to generate (near-equipartition) magnetic fields with coherence length 
L ~ c/ojp, and a field varying on such scale is expected to decay within a few skin-depths downstream. This 
suggests that the correlation length of the magnetic field far downstream must be much larger than the skin 
depth, L » c/u) p , perhaps even of the order of the distance from the shock. 

5. Observations of GRB afterglows imply that relativistic collisionless shocks produce a power- law distribution 
of high energy particles, d\ogn/d\oge = —2.2 ± 0.2. This is consistent with the results of diffusive (Fermi) 
acceleration of test particles, under the assumption of isotropic diffusion. However, a theory that self-consistently 
describes the acceleration of particles and the generation of electromagnetic waves is still missing. Since the 
high energy particles are likely to play an important role in the generation of the fields, a theory of collisionless 
shocks must provide a self-consistent description of particle acceleration, which depends on the scattering of 
these particles by magnetic fields, and field generation, which is likely driven by the accelerated particles. 

6. One of the important questions related to particle acceleration in GRBs is whether thes objects may produce 
the observed UHE, > 10 19 eV cosmic-rays. While there is evidence for an association between GRB and UHE 
cosmic-ray sources, the model for particle acceleration is incomplete. 

7. The search for a self-consistent theory of collisionless shocks has led to extensive numerical studies. 3D simula- 
tions provide compelling evidence that transverse electromagnetic instabilities can mediate relativistic collision- 
less shocks, that the upstream magnetic field is unimportant if the upstream is not highly magnetized, that the 
width of the shock transition region is a few xlO c/uo p , and that near-equipartition magnetic fields are generated 
in the downstream of relativistic shocks. 

8. These simulations are limited, however, to very small simulation boxes, and can reliably probe small length 
scales no larger than ~ 100 electron skin-depths, and short time scales no longer than ~ 100 electron plasma 
times. Obviously, the question of field survival and correlation length evolution on length scales 3> c/u p , which 
is the main challenge to our theoretical understanding, are not yet answered. Similarly, highly energetic particles 
cannot be contained in the small simulation boxes used, so Fermi-like acceleration processes are suppressed. 

9. The growth of the magnetic field correlation length to scales much larger than the skin depth suggests that 
downstream shock structure is self-similar. Assuming self-similarity allows one to draw important conclusion 
regarding the shock structure. In particular, the the magnetic field amplitude should scale with the distance 
D from the shock as B cx D 8B with — 1 < sb < 0, and the particle energy distribution should be a power-law, 
dn/de oc e~ p with p = 2/(sb + 1) (up to logarithmic corrections). If the low energy particles can be considered 
as infinitely conducting, the scaling is completely determined, sb — and p = 2. The self-similarity assumption 
may be tested by homogeneous plasma simulations. 
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